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ABSTRACT 

A series of three HST/STIS spectroscopic mappings, spaced approximately one year apart, reveal 
three partial arcs in [Fe II] and [Ni II] emissions moving outward from 77 Carinae. We identify these 
arcs with the shell-like structures, seen in the 3 D hydrodynamical simulations, formed by compression 
of the primary wind by the secondary wind during periastron passages. 

Subject headings: Stars: individual (77 Carinae) — stars: massive — binaries: general 


1. INTRODUCTION 

77 Carinae (hereafter 77 Car) is now recognized to 
be a massive binary system in a highly eccentric or- 
bit (e > 0 . 9 ). The primary star, 77 Car A, is in 
the luminous blue variable (LBV) phase, with a lu- 
minosity of 10 6,7 A© and mass of about 120 M© 
(Davidson & Humphreys 1997 ). It is surrounded by 
a massive stellar wind (M ~ 10 -3 M 0 yr _1 , v oo ~ 
420 kms -1 ; Groh et al. 2012) that interacts with a com- 
panion, 77 Car B, that has not been detected directly yet. 
The physical parameters of 77 Car B are, thus, inferred 
from X-ray (M ~ 10 -5 M 0 yr _1 , Vqq ~ 3000 kms -1 ; 
Pittard & Corcoran 2002) and photoionization modeling 
(~ 10 5 L 0 , 34000 K< T eff < 39000 K; Verner et al. 2005 ; 
Teodoro et al. 2008 ; Mehner et al. 2010 ). 

The orbital orientation is such that the orbital plane 
has an inclination of i ~ 140 °, with the semi- major 
axis pointing along position angle (P. A. 7 ) 314 ° and lon- 
gitude of periastron c 0 « 240 ° (Gull et al. 2009 , 2011 ; 
Madura et al. 2012). In this configuration, 77 Car B is 
in front of 77 Car A at apastron, and the orbital angular 
momentum vector is parallel to the symmetry axis of the 
bipolar Homunculus nebula that surrounds the system. 

Projected on the sky, at apastron, the secondary star 
directly illuminates and photo-ionizes the circumstellar 
material within 1200 AU of the binary system primar- 
ily to the northwest, where the material is primarily ap- 
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proaching the observer. The companion star spends most 
of the orbit near apastron, so that the undisturbed wind 
from 77 Car A is located on the far side of the system, 
red-shifted as seen from Earth, and largely unaffected by 
the ionizing radiation from 77 Car B. 

During periastron passage, 77 Car B approaches within 
1 — 2 AU of 77 Car A, leading to rapid changes in the 
shape of the wind-wind interaction surface and confine- 
ment of the far-UV flux from 77 Car B. Hydrodynamical 
simulations (Madura et al. 2012, T. Madura et al. 2013, 
in preparation) show that, after each periastron passage, 
a highly-distorted volume of hot, low density secondary 
wind pushes outward into the slow, high-density primary 
wind, leading to the formation of a thin, high density wall 
surrounding the lower density, trapped wind of 77 Car B. 
This dense wall is accelerated to velocities somewhat 
higher than the terminal velocity of 77 Car A and ex- 
pands both in the orbital plane and perpendicular to it, 
creating a thin, high-density sheet of trapped primary 
wind material. 

Using spectral imaging maps from the Space Tele- 
scope Imaging Spectrograph (S'T/S'), on board the Hub- 
ble Space Telescope ( HST ), we identify structures which 
are consistent with these walls of compressed primary 
wind material. We identify three partial arcs formed by 
the close passage of 77 Car B around the primary star 
over the last 3 orbital cycles, and using the STIS maps, 
we determine their space velocity and age. 

2. OBSERVATIONS, DATA REDUCTION, AND ANALYSIS 

Observations using HST/STIS were accomplished 
through three guest observer programs, 12013 (2010 Oc- 
tober 26), 12508 (2011 November 20), and 12750 (2012 
October 18), using the G430M and G750M gratings, cen- 
tered at 4706A and 7283A, respectively, in combination 
with the 52" x 1" aperture. An extended region, covering 
1" to 2" x 6.4", centered upon 77 Car, was mapped using 
0.05" spacing offsets at P. A. constrained by the space- 
craft. Exposures were approximately 30 seconds using 
CRSPLIT=2. An additional 7 position sub- map, also 
centered on 77 Car, employed shorter exposures to allow 
for potential saturation of continuum on the central po- 
sitions. 

Data were reduced with STIS GTO iDL-based 8 
tools similar to the standard pipeline reduction tools 

8 IDL is a trademark of Exelis Visual Information Solutions, Inc. 
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Fig. 1. — Iso- velocity images of [Fe II] A4815 for the observations 
of 2012 October 18. The velocities are indicated in each panel. The 
linear structures that appear in panels (g) and (h) are the result of 
a bad CCD column. The physical scale is based on an adopted dis- 
tance of 2350 pc (Smith 2006). We note that there are no extended 
structures at Doppler velocities higher than +440 kms -1 . 


(Bostroem & Proffitt 2011), but with improved ability 
to spatially align the spectra to sub-pixel accuracy. For 
each emission line of interest, continuum levels were fit- 
ted for each spatial position along the aperture and a re- 
sampled data cube was generated, with the coordinates 
of right ascension and declination at 0.05" spacing, and 
Doppler velocity at 20 kms -1 intervals. 

Since we are interested in structures formed in the 
primary wind around periastron, here we focus only on 
the following low-ionization lines: [Fe II] (Ao = 4729.39, 
4776.05, 4815.88, 7157.13) 9 and [Ni II] (A 0 = 7413.65). 
Observations from program 12013 (2010 October 26) did 
not include maps of [Fe II] A7157.13 and [Ni II] A7413.65. 

Visual examination of iso- velocity images in each data 
cube revealed multiple arcs in the low-ionization spectral 
lines used in this work. Sample images are presented in 
Figure 1 for the [Fe II] A4815 transition. These arcs 
are visible between P. A. 90° and 200°. We identified 
at least three sets of partial arcs, labelled Al, A2, and 
A3 in Figure 2. The two innermost arcs, Al and A2, 
are conspicuously detected in all of the transitions from 
[Fe II] and [Ni II]. The outermost arc, A3, is very faint. 

These arcs are not artifacts of the HST/STIS point 
spread function (PSF). The HST/STIS PSF diffrac- 
tive rings vary slowly with wavelength, are symmetric 
about the central core and have amplitudes substan- 
tially weaker than the partial rings apparent in the im- 
ages (Krist et al. 2011). As shown in the discussion be- 
low, these features correlate with multiple forbidden line 


9 All rest wavelengths, Ao, are in A, in vacuum (Zethson et al. 
2012 ). 
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Fig. 2. — Iso- velocity image of [Ni II] A7413 at +200 kms -1 , 
obtained in 2011 November 20. The arcs discussed in this work 
are indicated by the arrows and were labeled Al, A2, and A3. 
The dashed lines show the expected position of three spherical 
shells traveling at a constant speed of 475 kms -1 for 4.5, 10, and 
15.5 years, for the inner, middle, and outer arcs, respectively. To 
determine the expansion rate, we restricted our measurements to 
the eastern component of A2. 

structures seen at different wavelengths, expanding away 
from the central source. 

Comparison of the same iso- velocity image for different 
epochs (2010 October 26, 2011 November 20, and 2012 
October 18) revealed that Al, A2, and A3 move outward 
from the central source (Figure 3). Because A3 is faint 
and hard to see in all the images, and Al is affected by 
uncertainties in the continuum subtraction method close 
to the central source, we concentrate our analysis on A2. 
We note, however, that all three arcs seem to expand 
at the same rate, but some extended parts become very 
distorted or diluted, which makes the determination of 
accurate positions a difficult task. 

The position of the eastern portion of A2, relative to 
the central source, was measured using a radial flux pro- 
file cut along its direction (Figure 3). For each dataset 
(2010, 2011, and 2012), we extracted the flux profile be- 
tween 0 " < r < 0.9" and 100° < P. A. < 105°, where the 
emission from A2 is strongest. 

We used three different methods for measuring the 
peak position in the radial flux profile curve. The first is 
a simple gaussian model, which is a good approximation 
when the flux profile of the shell is somewhat symmet- 
ric. The second method relies on the barycenter of the 
flux profile, and the third, developed by Blais & Rioux 
(1986), is based on the analysis of the derivative of the 
observed profile. The position of the peak is the median 
value of these methods. 

A comprehensive study comparing several methods 
to determine peak positions with sub-pixel accuracy, 
including the three methods used here, was done by 
Fisher & Naidu (1996). The reader is referred to that pa- 
per for further details. Based upon the individual mea- 
surement errors, we estimate an accuracy of 0.25 pixel 
(0.013") for the position of the shell using the median 
value of the three methods mentioned before. The to- 
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Fig. 3. — Flux profiles and iso-velocity images of two forbidden lines: (a) and (e) are for [Ni II] A7413 and [Fe II] A4815, respectively, 
along the position indicated by the arrow in panel (6), which has a length of approximately 0.9 ,x and P. A. = 100°. The flux profiles 
were extracted from the iso- velocity images on the right. To illustrate the detected outward motion, in panels (c), (d), (/), (g), and (h) 
the curved solid line marks the position of A2 in 2010 October 26. For better visualization, the flux profiles were shifted in the vertical 
direction. 


tal error, however, may be larger than that because it 
is dominated by the dispersion in the measured position 
for various forbidden lines. 

3. RESULTS 

We note that A 2 is moving outwards with a proper 
motion /i = (0.05 ± 0.01)" yr -1 , on the plane of the sky 
(Figure 4(a)), a value suspiciously close to the pixel size. 
However, many other regions across the field of view are 
fixed in position from image to image, giving us con- 
fidence that the images are all spatially registered cor- 
rectly. Furthermore, the motion of portions of A 2 de- 
pends on the local Doppler velocity (cf. Figure 4(a)). 
Hence, the detected motion is not caused by a system- 
atic misalignment of the images. 

The projected radius (r p ), at a given Doppler velocity 
(I'd), of a shell moving with space velocity V, during a 
time £, is given by 


(475 ±29) kms -1 . The age of A 2 , t( A 2 ), in 2010, is thus 
obtained from < r p (A 2 ) > /V = (9.65 =b 0.55) yr, where 
< r p (A 2 ) >= 0.41" is the average projected distance of 
A 2 , relative to the central source, at that epoch. 

Assuming that A 1 is moving at the same constant 
space velocity as A 2 , then, for a given epoch, we have 


f p (A 1 ) = t( Al) 
r p (A 2 ) “ t( A2)’ 


(3) 


r p (t) = V t sin ^arccos 


Vb 

V 


(1) 


Taking the difference of the projected radius of A 2 , at 
two epochs separated by a time interval At = £2 — £ 1 , 
leads to 


r P {t 2 ) - r p (£ 1 ) = V At 


( ( V v 

( arccos ( — 




(2) 


Setting the observations made in 2010 Oct. 26 as base- 
line, we have At = 1.07 and 1.98 yr, for the two subse- 
quent datasets. Thus, we used the measured projected 
position for A 2 in the range —40 < Vd < +40 kms -1 
(this is the small angle regime, where Vd/V <C 1) 
to estimate the space velocity, which resulted in V = 


where r p (Al) is the average projected distance of Al 
from the central source, and £(Al) is its age. In 2010 
and 2011 , Al was too close to the central source to be 
accurately measured. Thus, we used the observations of 
2012 to measure the position and then estimate the age. 
In 2012 , we have r p (Al) « 0 . 22 ", resulting in £(Al) = 
(5.98 ± 0.92) yr. The relation described in equation (3) 
also applies to A3, resulting in £(A3) = (17.13 + 0.95) yr, 
since it was observed at r p (A3) ~ 0.6" in 2012 . The 
rather large errors are dominated by uncertainties in the 
position due to the proximity to the central source for 
Al, and the weak emission of A3. 

Correcting for the time interval between the 2012 and 
2010 observations, and comparing the ages in 2012 , we 
finally have £(A 1 ) = 5.98, £(A 2 ) = 11.63, and £(A3) = 
17.13 yr. This means that each shell is created approxi- 
mately 5.6 yr after the earlier one, showing that they are 
directly tied to the orbital period of the binary system 
(P = 5.54 yr; Damineli et al. 2008). 

To test for the consistency of our results across the 
range of Doppler velocities where A 2 is observed, we 
adopted a simplified model in which we consider A 2 to be 
part of a spherical shell centered on 77 Car, and allowed 
it to expand, at a constant space velocity of 475 kms -1 , 
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Fig. 4. — ( a ) Observed position of A2 as a function of Doppler velocity for the indicated epochs. The dashed lines show the projected 
position expected for a spherical shell expanding at 475 kms -1 during £l exp , £2 exp , and £3 exp . There is a clear departure from the spherical 
symmetry for red-shifted velocities > +100 kms -1 . (6) The dashed lines indicate the expected distance from the central source, as a 
function of time, of a spherical shell formed during periastron passage, and moving at 475 kms -1 . The word ’event’ corresponds to the 
periastron passage at the epoch indicated by the number in parenthesis. From this figure, Al was formed in 2009.03, A2 in 2003.49, and 
A3 in 1997.95. In 2010.8, there is no data for A3 because this arc can only be detected in the [Ni II] transition, and this transition was not 
included in the observations of that year. 


during the corresponding age of A2 at each epoch of ob- 
servations (tlexpj t2 exp , and £3 ex p in Figure 4(a)). 

At line of sight velocities near zero, a spherical model 
reproduces the observations very well, with standard 
deviations < 0.015". However, it fails to match the 
observed position of A2 at red-shifted velocities > 
+100 kms -1 ; A2 moves slower than predicted. The final 
standard deviation for the spherical model is 0.05", in 
the range from —40 to +220 kms -1 . 

Figure 4(b) shows the comparison between the ob- 
served and the expected motion of shells formed during 
various periastron passages, and moving at 475 kms -1 . 
Despite showing an offset between the expected and ob- 
served position at a given time, their slope are the same, 
which gives support to our derived space velocity and 
the assumption that the shells are expanding at the same 
rate. 

4. DISCUSSION AND CONCLUSIONS 

We used the critical density for each forbidden line, cal- 
culated using the CHIANTI atomic database (Dere et al. 
1997; Landi et al. 2013), to infer the typical local den- 
sities of the arcs. The outermost arc, A3, seen in the 
[Ni II] transitions but weakly in the [Fe II], must have a 
density in the range 0.8 — 4 x 10 6 cm -3 . The strong [Fe II] 
emission from Al and A2, indicates that their densities 
must be > 4 x 10 6 cm -3 . 

Assuming that A2 is half a spherical shell ( cf . Figure 5) 
with average radius of 0.4" (projected distance from the 
central source, at zero Doppler velocity, in 2010) and 
thickness < 0.05" (since we cannot spatially resolve the 
arcs), and using the critical density calculated before, 
we estimate a total mass of 2.6 x 10 -3 M 0 , for solar 
metallicity. 

Models for the primary star predict that its wind den- 
sity should drop below 10 6 cm -3 for projected distances 
greater than 0.2" (Hillier et al. 2001; Groh et al. 2012). 
Since we detected arcs exceeding 10 6 cm -3 out to 0.6", 


this is naturally explained by the primary wind being 
compressed - a factor of 30 - by the secondary wind 
during periastron passage. This is also supported by the 
fact that the average age difference between the arcs is 
approximately 5.6 yr, nicely matching the binary sys- 
tem’s period. 

The primary mass-loss rate required to produce a shell 
with the same mass and geometry as A2, over a time in- 
terval equal to the binary’s period, is 9.3 x 10 -4 M 0 yr -1 , 
in agreement with the results from radiative transfer 
spectral modeling (Hillier et al. 2001, 2006; Groh et al. 
2012). We note, however, that our result must not be 
taken as a strict lower or upper limit, since the actual ge- 
ometry and thickness of A2 may be different from what 
we assumed. 

Simulations using 3D smoothed particle hydrodynam- 
ics (3D SPH) reveal compressed regions of the primary 
wind formed during the periastron passages and main- 
tained for multiple cycles by the hot, low density, resid- 
ual secondary wind (Figure 5). These simulations pre- 
dict that these compressed regions accelerate from the 
primary wind terminal velocity (uoo = 420 kms -1 ; 
Groh et al. 2012) to about 460 — 470 kms -1 (T. Madura 
et al. 2013, in prep.), which is consistent with our derived 
expansion velocity of 475 kms -1 for the shells. 

The departure from spherical symmetry for high red- 
shifted Doppler velocities suggests that the shells do not 
expand with the same velocity in all directions. This 
is explained by the fact that 3D SPH simulations show 
that material located in the orbital plane moves faster 
than material located off the plane. Thus, the observed 
Doppler velocity is a combination of velocities from dif- 
ferent emitting regions along a specific viewing angle, 
which, in turn, depends upon the orientation of the shells 
on the sky, ultimately defined by the orbital orientation. 

During periastron passages, SPH 3D simulations show 
that the secondary star creates a cavity deep within the 
primary wind, pushing a shell outwards. Since the low 
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Fig. 5. — 3D SPH simulations showing the formation and geometry of high density regions due to the compression of the primary wind 
by the secondary wind, in two different regimes for the primary mass-loss rate: 8.5 X 10 -4 Mq yr -1 (left) and 2.5 x 10 -4 Mq yr -1 (right). 
Upper panels : slices in the orbital plane (xy). Lower panels : slices perpendicular to the orbital plane (yz). In the high mass-loss rate 
regime, there is a gap in the polar region that is not observed in the lower mass-loss rate regime. The lower and upper limits of the density 
scale correspond to the critical density values for [Ni II] and [Fe II], respectively. 


density secondary wind is moving at « 3000 kms -1 , the 
interaction with the high density, relatively slow primary 
wind results in the formation of a high density, thin shell 
accelerated to velocities slightly higher than the terminal 
velocity of the primary wind. Hence, the geometry of the 
shell is strongly affected by the physical parameters of 
the binary system, and the projected appearance on the 
sky is defined by the inclination of the orbital plane. 

The latitudinal extent of the shells is controlled by the 
asymptotic opening angle, 0 oo, of the wind- wind collision 
region, given by (Usov 1992) 


0 


oo 


27 r 

Y 



(4) 


where 77 = Mb ^ooB / Ma f 00 A is the ratio of the sec- 
ondary to the primary wind momentum, M and are 
the mass- loss rate and terminal velocity of the wind, and 
the subscript A and B refer to the primary and secondary 
star, respectively. 


Assuming fixed parameters for the secondary wind, a 
large primary mass loss rate decreases the bow shock 
opening angle #oo, so that the compressed primary wind 
material is trapped near the orbital plane (Figure 5(c)); 
a low mass-loss rate increases #oo, and the compressed 
primary wind can extend well off the orbital plane (Fig- 
ure 5(d)). This means that the geometry of the circum- 
stellar arcs is an excellent diagnostic of the primary mass- 
loss rate: a lower mass loss rates produces a nearly con- 
tinuous arc seen in projection around the star, while a 
higher mass loss rate produces a shorter, broken arc. 

The existence of a discontinuity in the arcs that we 
observe, more pronounced in Al as a lack of emission at 
P. A.« 135°, yields a lower limit to the mass-loss rate of 
the primary star. SPH simulations in 3D show that such 
broken rings are only apparent if the primary mass-loss 
rate is above 2.5 x 10 -4 M© yr -1 , which suggests that, at 
least for the last event, the primary mass- loss rate must 
have exceeded this value. 

Our results also constrain the orbital orientation, with 
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periastron occuring between 200° < uj < 270°, z.e., dur- 
ing periastron, the primary star is between the observer 
and the secondary star. The opposite configuration, 
where the secondary, at periastron, is between the ob- 
server and the primary star, would produce blue-shifted 
shells seen to the northwest of the central source, rather 
than the red-shifted features to the southeast that we 
observe. 

Finally, we note that our results assume a constant 
shell space velocity. However, as can be seen in Fig- 
ure 4(b), there is an offset of 0.07" between the expected 
and the observed distance from the central source as a 
function of time. We know, from the comparison between 
the derived space velocity (475 kms -1 ) and the terminal 
velocity of the primary wind (420 kms -1 ), that this off- 
set is due to the acceleration that these shells experienced 
at some point during the early stages of their formation. 
Thus, the ages derived in this letter are overestimated, 
i.e the shells are younger. Nevertheless, the difference of 
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5.6 yr between them remains the same, validating their 
connection with the orbital period. 
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